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ABSTRACT 

We present the X-ray spectral results from the longest X-ray multi-mirror mission- 
Newton observation, 133 ks, of the low luminosity active galactic nucleus NGC 7213. 
The hardness ratio analysis of the X-ray light curves discloses a rather constant X- 
ray spectral shape, at least for the observed exposure time, enabling us to perform 
X-ray spectral studies using the total observed spectrum. Apart from a neutral Fe 
Ka emission line, we also detect narrow emission lines from the ionised iron species, 
Fe XXV and Fe xxvi. Our analysis suggests that the neutral Fe Ka originates from 
a Compton-thin reflector, while the gas responsible for the high ionisation lines is 
collisionally excited. The overall spectrum, in the 0.3-10 keV energy band, registered 
by the European Photon Imaging Camera, can be modelled by a power-law component 
(with a slope of F ~ 1.9) plus two thermal components at 0.36 and 8.84 keV. The low- 
energy thermal component is entirely consistent with the X-ray spectral data obtained 
by the Reflection Grating Spectrometer between 0.35-1.8 keV. 
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1 INTRODUCTION 

NGC 7213 (z=0. 005839) is a face on Sa galaxy hosting an 
active galactic n ucleus (AGN). It has been classified as a 
type I Seyfert bv lPhillipsI { 19791 ). based on its linewidth 
(full width at zero intensity of 13000 km s~^), but also as 
a low-ionisation nuclear em ission line region (LINER) by 

iFilippenko fc HalpernI (|l984l '). based on the study of a vari- 
ety of optical emission lines which were observed to have a 
' full width at half maximum of 200 to 2000 km s~^. 

NGC 7213 hosts a blac k hole with a mass of about 
10** M0 (|Woo fc Urrvll2002l '). It h as a bolometric luminos- 
ity o f Lboi = 1.7 X 10*^ ergs s ~^ l|Emmanoulopoulos et al.l 

|2012| ). yielding a rather low accretion rate of 0.0014 (i.e. 
0.14 per cent) of the Eddington limit being intermedi- 
ate between those usually found in type I Seyfert galax- 
ies and liner's. NGC 7213 exhibits also intermediate ra- 
dio properties between those of radio-loud and radio-quiet 
AGN. Its long-term X-ray spectral shape evolution, stud- 



* Based on observations obtained with XMM-Newton, an ESA 
science mission with instruments and contributions directly 
funded by ESA Member States and NASA, 
t E-mail: D.Emmanoulopoulos@soton.ac.uk 



ied with the use of Rossi X-ray Timtng Explorer (RXTE) 
data, show a clear 'harde r when brighter' X-ray behaviour 
jEmmanoulopoulos et all |2012| ). Both these facts make 
NGC 7213 an extragalactic analogue of the Galactic 'hard 
state' sources. At the same time its radio and X-ray lumi- 
nosities can be fitted with the global fundamental plane be- 
ing als o consistent wit h a 'hard state' identification of this 
AGN (|Bell et al.ll2oTlh . 

NGC 7 213 has been o bserv ed several times in the X- 
ray band. iBianchi et al.l (|2003l ) used simultaneously ob- 
tained observations from the European Photon Imaging 
Camera (EPIC), on board X-ray multi-mirror mission- 
Newton (XMM-Newton) observatory, and BeppoSAX satel- 
lites. They found no evidence of a Compton reflection com- 
ponent and detected an 'iron line complex' between the 
energies of 6.5 and 7 keV, which could be explained in 
terms of three narrow emission lines, Fe Ka, Fe xxv and 
Fe xxvi, produced either in the Compton-thin torus or in 
the broad line region (BLR). Based on the Reflection Grat- 
ing Spectrometer (RGS) data, of the same X MM-Newton 
pointing (around 46 ks). IStarling et all |200^ did not find 
any absorption features (like those typically observed in 
type I Seyfert galaxies). Instead, they detected emission 
lines which could arise from a collisionally ionised ther- 
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mal plasma with a temperature of about 0.18 keV. Fur- 
ther observations using the high-energy tran smission grat- 
ing ( HETG), on board Chandra observatory ijBianchi et ahl 
|2008| ). suggested that the neutral Fe Ka emission line can 
be produced in the BLR since its measured full width at 
half-maximum of 2400lg|!,Q° km s~^ was fully consistent 
with that of the line (2640lgo° km s"^) in this ob- 
ject. The same authors found out that the dominant com- 
ponent of the Fe xxv triplet is the res onant line suggest- 
ing an origin in coUisionally ion ised gas (|Porauet fc Dubaul 
l2000l : lBautista fc KallmanH 2000'') . A broad-band X-ray study 
of NGC 7213, 0. 6-150 keV using ob servations from Suzaku 
and Swtft BAT l|Lobban et al.ll2010l ), derived similar results 
with respect to the origin of the neutral Fe Ka line and 
the absence of the Compton reflection component suggest- 
ing tha t the inner, opt i cally thick accretion disc, is probably 
absent. iLobban et al.l (|2010l ) also found that the dominant 
component of the Fe xxv is consistent with that of forbid- 
den transition, originating from highly ionised pl asma and 
not fr om coUisionally ionised gas as found by Bia nchi et al.l 
ll2008D . 

In this paper we report the analysis results of the longest 
XMM-Newton observation, around 133 ks, of NGC 7213 per- 
formed in November 2009. Initially, in Sect.[2]we present the 
data-reduction procedures for the EPIC X-ray data, consist- 
ing of the pn-charge coupled device (pn-CCD) and the two 
metal oxide semi-conductor (MOS)-CCDs. In the same sec- 
tion, we present the data-reduction details for the soft X- 
ray data from the RGS. Then, in Sect. [3] we show the EPIC 
light curve products together with a hardness-ratio analy- 
sis. In Sect.m we perform the spectral analysis of the EPIC 
data and in Sect. [5] we investigate the spectral properties 
of the RGS data. A discussion together with a summary of 
our results can be found in Sect. [S] The cosmological pa- 
rameters used throughout this paper are: Ho = 70 km 
Mpc~^, Qa ~ 0.73 and fim — 0.27, yielding a luminosity dis- 
tance to NGC 7213 of 22.12 Mpc (for a corrected redshift, 
«corr.3K = 0.005145 to the reference frame defined by the 3 
K cosmic microwave background radiation). This value of 
the luminosity distance appears to be fully co nsistent wit h 
the redshift-independent TuUy-Fisher distance (|Tu11vII 19881 ). 
Throughout this paper the error bars in the figures corre- 
spond to the 68 per cent confidence limits of the correspond- 
ing plot points. 



2 OBSERVATIONS AND DATA-REDUCTION 

2.1 EPIC-MOS and EPIC-pn data-reduction 

NGC 7213 was observed by XMM-Newton (Obs-ID: 
0605800301) from 2009 November 11, 21:05:46 (UTC), to 
2009 November 13, 09:54:25 (UTC) (on-time: 132519 s). The 
pn camera was operated in Prime Small Window mode and 
the two MOS cameras in Prime Partial W2 (i.e. Small Win- 
dow) mode. Medium-thickness aluminised optical blocking 
filters were used for all EPIC cameras to reduce the contam- 
ination of the X-rays from infrared, visible, and ultra-violet 
light. 

The EPIC raw-data are reduced wi th the XMM-Newton 
SCIENTIFIC ANALYSIS SYSTEM (SAS) (jCabriel et al.l I2OO4I ') 
version 12.0.1. After reprocessing the pn and the two MOS 



data-sets with the epcham and the emchatn SAS-tools re- 
spectively, we perform a thorough check for pile-up using 
the task epatplot. Neither the pn nor the two MOS data 
sets appear to suffer from pile-up, therefore count-rates for 
the source and the background are extracted from circu- 
lar regions of radii (in arcsec): 37.75, 45 (pn), 39.5, 179 
(MOSl) and 43.5, 200 (MOS 2). The corresponding centres 
of the circular regions for the source and the background 
are the following (in camera coordinates: [RAWX,RAWY]): 
pn -[27660,27023], [28204,23151] 3.3 arcmin of the source 
to the south east, (MOS 1) -[27663,27029], [15947,27089] 9.6 
arcmin of the source to the west-, (MOS 2) -[27653,27031], 
[18569,18387] 10.7 arcmin of the source to the south west-. 
Note that these apertures yield the optimum signal-to-noise 
ratio above 7 keV. We verify that the resulting light curves 
are not affected by pile-up problems. 

For the extraction of the light curves at a given en- 
ergy range, we select events that are detected up to quadru- 
ple pixel-pattern on the CCDs i.e. PATTERN<12. The cor- 
rected background-subtracted light curves of the source are 
produced using the SAS-tool epiclccorr. Based on the pn 
background light curve, in the 0.3-10 keV energy range 
(Fig. m grey-points corresponding to the right-axis), an in- 
creased background activity was registered during the first 
7 ks and the last 8.9 ks of the observations and thus, data 
from those periods are disregarded from our analysi^J. This 
leaves us with 116619 s clean exposure-time. 

For the production of the spectra we use the SAS-task 
evselect and we select for the pn camera the pulse-invariant 
(PI) channels from to 20479, having a spectral bin size of 
5 eV. For the two MOS cameras we select the PI channels 
from to 11999 with a spectral bin size of 15 eV. For the 
spectral analysis we set more stringent filtering criteria from 
those in the case of the light curves. For the pn data we allow 
events that are detected up to double pixel-pattern on the 
CCD i.e. PATTERN<4 and we exclude all the events that 
are at the edge of a CCD and at the edge to a bad pixel 
i.e. FLAG=0. For the MOS data we allow, as in the case 
of the light curves, PATTERN<12 but we restrict to those 
events having also FLAG=0. After filtering the events, we 
calculate the area of the source and the background circular 
regions using the task backscale and finally we compute the 
detectors' response matrices and the effective areas using 
rmfgen and arfgen, respectively. 

The X-ray spect ral fitting anal ysis has been performed 
by the xspec task l|Arnaudl [l996l ') version 12.7.1, of the 
HEASOFT version 6.12. In order to have similar signal-to- 
noise ratios among different spectral bins we group the pulse- 
invariant (PI) channels of the source-spe ctra, using the task 
grppha of the ftools (|Blackburnlll995l ') version 6.10. The 
MOS source-spectra are grouped in bins of four (Ai5=60 
eV) and eight (A_E=120 eV) spectral channels for the en- 
ergy ranges of 0.3-8 keV and 8-10 keV, respectively. The 
pn source-spectra are grouped in bins of four {AE =20 eV) 



Despite the fact that the 'high' background-states are still min- 
imal, in terms of count-rate, with respect to the source count-rate 
in 0.3-10 keV (Fig. [T] black-points corresponding to the left- 
axis), they consist mainly of events having energies more than 
2 keV. At these energies the source count-rate is of the order of 
1.5 counts-s~-^ . 
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Figure 1. The NGC 7213 EPIC pn and MOS combined light 
curve (black-points, left-axis) and the background pn light curve 
(grey-points, right-axis) in the observed 0.3-10 keV energy-band. 
The background shows increased activity during the first 7 ks and 
the last 8.9 ks. 



spectral channels for the energy range of 0.3-10 keV. All the 
resulting grouped spectral bins contain at least 20 photons. 



2.2 RGS Data Reduction 

The RGS raw-data are processed by following the standard 
data-reduction threads from SAS. We use the tool rgsproc 
to extract calibrated first order spectra and responses for 
the RGS 1 and RGS 2 cameras. Just like the EPIC cameras 
the RGS observations can be afi'ected by high particle back- 
ground periods during parts of the XMM-Newton orbits, 
mostly caused by Solar activity. The high energy band is 
the one most affected by background flares, and high energy 
RGS photons are dispersed over the CCD-9 chip. Given this 
fact, we extract the background light curve of the CCD-9 
chip and select, as good-time-intervals of the processed ob- 
servations, only those during which the background count 
rate deviates by less than 2 standard deviations from the 
average background count rate of each observation. The fi- 
nal (i.e. after cleaning for high-background time intervals) 
exposures of the RGS 1 and RGS 2 spectra are 125.3 and 
125.6 ks, respectively. 

Both RGS 1 and RGS 2 spectra are grouped at a resolu- 
tion of 160 mA for visual purposes (i.e. for the figures where 
we present the RGS spectra). This choice corresponds to 
around 3.2 eV at 0.5 keV, and is equivalent to sub-sampling 
the effective RGS resolution by a factor of 2. Since the signal- 
to-noise ratio per resolution element is too low in the con- 
tinuum, the chosen grouping helps to pinpoint the various 
emission and absorption features easier through optical in- 
spection. For model fitting purposes, we bin the RGS spectra 
at a resolution of 0.015 A, i.e. we oversample the effective 
RGS resolution by a factor of 5. 

Due to failures of two difi^erent read-out detector chips, 
early in the life of the mission, both the RGS 1 and the RGS 2 
lack response in two different spectral intervals of 0.9-1.2 
keV and 0.5-0.7 keV, respectively. We therefore consider the 
following two spectral intervals for spectral fitting purposes: 
[(0.35 - 0.9) U (1.2 - 1.8)] keV and [(0.35 - 0.5) U (0.7 - 1.8)] 
keV, for RGS 1 and RGS 2, respectively. The upper energy 



^ limit of 1.8 KeV is based on the fact that above this energy 
o the RGS effective area drops from 70-100 cm^ to 30 cm^ and 
iq as a result the observed spectra become very noisy (carrying 
§ no extra value for model fitting purposes). 
2. The net counts per bin (i.e. after background subtrac- 
g tion) at the corresponding spectral intervals are between 100 
S and 300, yielding a signal to noise ratio between 10 and 17 
" per bin for RGS 1 and RGS 2, respectively. Outside this band 
the photon count-rate of the background spectrum becomes 
rT comparable to the background-subtracted source count-rate, 
c For the spectral fi tting procedures we have used the 
^ Sherpa fitting package (|Freeman at al.|[200ll ). of the chan- 

DRA INTERACTIVE ANALYSIS OF OBSERVATIONS (CIAO) Soft- 

ware version 4.4, to fit simultaneously the RGS 1 and RGS 2 
spectra of NGC 7213. 



3 X-RAY LIGHT CURVES 
3.1 XMM-Newton light curves 

In order to produce the combined X-ray light curves from 
the EPIC-pn and the two EPIC-MOS cameras we estimate, 
in 100 s bins, the total count rate and the one standard de- 
viation error bar, as the Euclidean norm of individual errors 
(i.e. the square root of the sum of the squared errors) from 
all three detectors, using the ftool Icmath. The combined 
final X-ray light curves in 0.5-1.5 keV, 2.5-4 keV, and 5-10 
keV energy ranges, are shown in Fig. (2] A very small steady 
count-rate increase can be noticed in all three energy-bands, 
of the order of 6.7 ± 0.5, 5.7 ± 1.3 and 6.4 ± 1.7 per cent, 
respectively. Additionally, the fractional va riability ampli- 
tudes (corrected for the photon noise e.g. IVaughan et al.l 
l2003h . for the three bands are rather small, 1.7±0.3, 2.5±1.6 
and less than 0.8 per cen10, respectively. 



3.2 XMM-Newton hardness-ratio analysis 

After binning the light curves in bins of 5 ks, we estimate 
the hardness-ratios (5-10 keV)/(0.5-1.5 keV) (HRl) and 
(5-10 keV)/(2.5-4 keV) (HR2) versus the overall count- 
rate in 0.5-10 keV (Fig. [3]). There is no direct evidence 
of an increasing or decreasing trend which would imply 
that the X-ray spectrum of the source becomes softer or 
harder when the source gets brighter. The error bars cor- 
respond to one standard deviation and they are derived 
following the standard erro r propagation procedure (e.g. 
iBevington fc Robinsonlll992h . 

In order to check statistically HR 1 for a potential trend, 
we fit to the 22 hardness-ratio points a linear model and we 
check how significantly the value of the slope differs from 
zero. For HR 1 the best-fitting model gives a slope of ( — 1.8± 
4.3) X 10"^ with a 95 per cent confidence interval (c.i.) of 
(-1.1 X 10"^7.2 X 10"^lf|. Since the value of f-statistic that 
we get from the data-set is 3.37, the probability of getting 
such a value by chance alone is 3.05 x 10"'^. Similarly, for 
HR2 the best-fitting model gives a slope of (1.2±3.2) x 10"^ 
having a 95 per cent c.i. of (-5.5 x 10"^, 8.0 x 10"^). The 



^ The excess variance of the light curve in the 5-10 keV energy 
band is negative and thus we quote the 90 per cent upper limit. 
The value of t-statistic is t20, 0.025 = 2.09. 
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Figure 2. The combined EPIC-pn,-MOS 1,-MOS 2 X-ray light 
curves of NGC 7213 in the 0.5-1.5 keV, 2.5-4 keV, and 5-10 keV 
energy-bands, in bins of 100 s. Note the difference in the count- 
rate scale on the vertical axis. 

value of t-statistic for this data-set is 2.36 yielding a chance 
coincidence probability of 2.8 x 10~^. 

These results imply that NGC 7213 does not show sig- 
nificant X-ray spectral variations, as a function of time, and 
thus we can study the average overall X-ray spectra without 
being affected from temporal spectral variations. 



4 EPIC X-RAY SPECTRAL ANALYSIS 

During the X-ray fitting procedures of the EPIC data we 
keep all the model-component parameters between the two 
MOS and the pn spectra tied together except from the nor- 
malisations (for the two MOS spectra they are still tied to- 
gether), which are always left varying freely, unless otherwise 
stated. The model-flux is modifled using the photoelectric 
absorption model wabs having a fixed interstellar column- 
density equal to the measured hydrogen Galactic column of 
the neutral interstellar medium along this line of sight of 
A^H = 1.1 X 10^° cm~ ^ (estimated using the ftool nH, after 
iKalberla et"alll2005l ). 

Note, that the actual energies of the emission lines 
as well as the best-fit values for their centroid energies 
are given in the source's rest-frame. The horizontal axis 
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Figure 3. Hardness-ratios in (5-10 keV)/(0.5-1.5 keV) and (5- 
10 keV)/(2.5-4 keV) versus the total count-rate in 0.5-10 keV, 
binned in 5 ks time-bins. The dotted-line represents the best- 
fitting linear model. 



in the spectral plots always refers to the energies in the 
observer's frame. For the X-ray spectral results, we use 
the fit statistic, with data variance, tracing the best- 
fit mo del parameters using the Lev enberg-Marquardt algo- 
rithm (|Bevington fc Robinson|[l992h . The errors on the lat- 
ter indicate their 90 per cent confidence range, correspond- 
ing to a Ax^ of 2.706, unless otherwise stated. Finally, the 
best-fit parameters for the following X-ray spectral models, 
at 3-10 and 0.3-10 keV, are given in Table [T] 



4.1 A first look at the broad-band X-ray spectrum 

In order to get a phenomenological view of the overall EPIC 
X-ray spectra we initially fit a redshifted power-law (Model 
1) to the data in the 3-5 keV and 7-10 keV energy-bands. 
The x^ value of the best-fit model is equal to 404.57 for 387 
d.o.f. (0.259) null-hypothesis probability (NHP)) yielding a 
best-fit photon spectral index of F = 1.66 ± 0.02. Figure |4] 
shows this best-fit model extrapolated over the energy range 
of 0.3-10 keV together with the respective ratio plot. 

From the ratio plot of Fig. [4] we can readily distinguish 
three very distinct spectral features. Firstly, below 2 keV, 
we notice a fiux excess on top of the best-fitting power- 
law extrapolation. Then, we can discern a strong emission 
line appearing around 6.36 keV and finally several excess 
features between 6.4 and 7.2 keV. 

For clarity reasons, in Fig. [5] we show a blow-up of the 
same ratio plot, for the region 5.9-7.2 keV. We plot only 
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Figure 4. Model 1: Redshifted power-law fit in the observed 3- 
5 keV and 7-10 keV energy ranges, indicated by the grey-areas. 
The best-fit model is extrapolated over the energy range of 0.3—10 
keV. 
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Figure 5. The EPlC-pn spectral data between 5.9 and 7.2 keV 
together with the expected positions of some Fe species indicated 
by the vertical dashed lines. 



the ratio points of the pn data, together with the position 
of the following iron emission lines that are frequently ob- 
served/assumed in the AGN X-ray spectra: a neutral Fe 
Kq (6.4 keV), He-like Fe xxv- [forbidden (f), intercombina- 
tion (i) and resonance (r)] (6.637, 6.668, and 6.700 keV), H- 
like Fe xxvi [ls-2p] (6.966 keV) and a neutral Fe K/3 (7.06 
keV). 



4.2 Spectral fitting in the 3—10 keV energy range 

Since Model 1 (Fig. U) gives an adequate description of the 
data in the 3-5 keV and 7-10 keV energy ranges, we keep it 
as a baseline model, for the following spectral analysis over 
the 3-10 keV energy range. As we discuss below, the ex- 
cess emission between 6.4 and 7.2 keV is modelled as an en- 
semble of individual narrow Gaussian function having fixed 



the standard deviation at 22 eV. This value corresponds 
t o a full width half maximum (FW HM) of 2400 km s"^ 
(iBianchi et all l2008l : IShu et all bOlOl ) and it is well below 
the actual energy resolution of the EPIC instruments (i.e. 
around 150 eV FWHM at 6.4 keV for both MOS and pn). 



4-2.1 Iron line species 

Initially, in order to account for the observed emission 
line at 6.37 keV, which appears to correspond to the Fe 
Ka line from neutral iron (having a centroid energy of 
6.4 keV), we consider a redshifted narrow Gaussian func- 
tion. For consistency reasons with the theoretical predic- 
tions (jKaastra & Mewe 1993), the presence of a neutral Fe 
Ka requires the existence of an associated Fe K/3 line. There- 
fore, we add another redshifted narrow Gaussian function 
to account for the Fe K/3 line, fixing its centroid energy at 
7.06 keV and its normalisation at 14 per cent of tha t of the 
Fe Ka 's normal isation value l|Molendi et al.ll200j ). Using 
Chandra-HETG, iBianchi et al. I (|2009l ') found no evidence of 
a Fe K/3 line yielding a 90 per cent upper confidence limit for 
the for its normalisation of 18 per cent of the Fe Ka 's nor- 
mahzation, being consistent with the theoretically expected 
value. However, the presence of a narrow K/3 line is not sta- 
tistically required from our spectral data (i.e. its addition to 
the models described below does not significantly improve 
the quality of the model fits). However, for physical reasons, 
given the presence of the strong neutral iron Ka emission 
line in the spectrum of the source, we decide to keep this 
emission line in our X-ray spectral fits (note that we veri- 
fied that our model fit results are not significantly altered if 
we exclude this component). 

The best-fitting model (Model 2, Fig. |S1 top panel) has 
a value of 696.13 for 547 d.o.f. (1.5 x 10"^ NHP) and its 
best-fitting values for the Fe Ka's centroid energy and the 
power-law's photon index are 6.42±0.01 keV and 1.63±0.02, 
respectively. From the corresponding ratio plot (Fig. [S] top 
panel) we can readily see that the ratios around 6.37 keV 
and 7.02 keV are consistent with unity but there is still a 
great deal of spectral excess between the observed energy 
range of 6.5-7 keV. 

In order to model this excess we consider additionally 
a third redshifted narrow Gaussian function (Model 3). The 
best-fitting Model 3 gives a considerably better value of 
617.05 for 544 d.o.f. (1.6 x 10"^ NHP) yielding for the new 
component a best-fitting centroid energy of 6.731q q5 keV. 
The centroid energy of the Fe Ka line as well as the photon 
index of the power-law are still consistent with the previ- 
ously derived values (i.e. Model 2). The bottom panel of 
Fig. |6] shows the best-fitting model together with the cor- 
responding ratio plot, the latter, showing an excess around 
6.8 keV. 

In order to model this excess we add a fourth redshifted 
narrow Gaussian function (Model 4, Fig. [Tjl . The addition of 
this narrow line component ameliorates significantly the 
value of the best-fitting Model 4, 593.41 for 541 d.o.f. (0.06 
NHP), yielding a best-fitting centroid energy of 6.961q'q5 
keV. The best-fitting centroid energy of the third narrow 
emission line (i.e. the one introduced in Model 3) is now 
6.71 ± 0.02 keV and the best-fitting values for the Fe Ka's 
centroid energy and power-law photon index remain exactly 
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Figure 6. Line spectral fits in the observed 3-10 keV energy 
range. The vertical labels indicate the energy in the source's 
rest-frame. (Top panel) Model 2: Two Fe lines [Fe K« (free), 
Fe K/3 (fixed)] ;a redshifted power-law. (Bottom panel) Model 3: 
Three Fe lines [Fe Ka (free), Fe xxv-r (free), Fe K/3 (fixed)] ;a 
redshifted power-law. 



the same with those ones derived previously from the best- 
fitting Model 3. 

The narrow emission line, having a best-fitting centroid 
energy of 6.96 keV, should correspond to the ls-2p doublet 
from H-like Fe xxvi. The other emission line, having a best- 
fitting centroid energy of 6.71 keV, is consistent with the He- 
hke Fe xxv-r since its 99 per cent c.i. is (6.667-6.741) keV 
(Fig. [H]) . Using the best-fitting results of Model 4, we can now 
estimate the equivalent widths (EWs) for the various iron 
line species, having as an underlying continuum the power- 
law with photon index 1.66. For the Fe Ka and the Fe K/3 we 
get 126^27 and 23^5 eV respectively, and for the Fe xxvi 
and Fe xxv-r we get 23I23 ^^'^ ^Sljo eV respectively. 
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Figure 7. Line spectral fits in the observed 3-10 keV energy 
range for Model 4. The vertical labels indicate the energy in the 
source's rest-frame. Four Fe lines [Fe Ka (free), Fe xxv-r (free), 
Fe xxvi (free), Fe K/3 (fixed)] ;a redshifted power-law. 



4.2.2 Origin of Fe Ka line 

Next we consider whether the existence of a distant 
Compton-thick refiector can further improve the fit pro- 
viding at the same time a natural explanation for the Fe 
Ka line. To test that we fix all the iron line components to 
their best-fitting values from Model 4, apart from that of Fe 
Kq line, and we include to that the additive xspec model 
pexrav (jMaedziarz fc Zdziarski|[l995h . During the fit, we tie 
the photon index of the pexrav-model to that of the power- 
law continuum, we fix the abundance of elements heavier 
than helium to that of the solar value and we fix the cut-off 
energy to 1000 keV. Then, we consider an ensemble of in- 
clinations 30°, 50° and 80° as well as iron abundances 0.6, 
0.8 and 1 and we estimate the best-fitting reflection scaling 
factor R for each case. For all the cases, the best-fit pexrav 
normalisation is equal to zero having a value of 589.65 for 
541 d.o.f. (0.07 NHP). The reduced x^s between this model 
and Model 4 (1.09 and 1.10 respectively) imply that such a 
Compton-thick distant refiection component is not required 
from our spectral data. 

Note that by excluding all the iron line species (i.e. 
Gaussian components) and fit the spectral data, with the 
above-mentioned set of pexrav models, we get as a best fit 
model the one having a frozen inclination of 50° and an iron 
abundance of 0.8. This model yields a value of 1009.11 for 
548 d.o.f. (0 NHP) being worse than that of Model 4 and thus 
making the presence of the iron line components statistically 
more necessary than that of a Compton-thick distant re- 
fiection component. Therefore, in the following broad-band 
spectral fits we exclude the pexrav-model from our analysis. 
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Figure 8. The versus the centroid rest energy of the 6.71 keV 
emission line of Model 4. The thin vertical grey lines indicate the 
rest energies of the Fe XXV triplet components. The horizontal 
dot-dashed lines indicate the 90, 95 and 99 per cent confidence 
ranges. 



4.3 The 0.3-10 keV EPIC-pn spectrum: Thermal 
plasma components 

As we have seen from the ratio plot of Fig. U] NGC 7213 
exhibits below 2 keV an excess flux (i.e. additional to the 
power-law component) known as 'soft excess' emission. 

In order to fit the 'soft excess' we extend the best- 
fit Model 5 down to 0.3 keV a nd we incliide the additive 
XSP EC thermal model mekal (,Mewe et al.ll 19851 : iRavmondl 
I2OO9I ) incorporating the emission spectra from hot diffuse gas 
(thermal plasma). During the fit we fix the energies of the 
iron line species to their best fitting values as derived from 
Model 4 and for the mekal-model we freeze both the Hydro- 
gen volume density and the chemical abundance of metals 
to unity and to that of Solar abundance, respectively. The 
best-fitting model (Model 5a, Fig. (5] top panel) has a 
of 896.44 for 764 d.o.f. (6.26 x 10"" NHP) yielding a best- 
fitting plasma temperature of ksT = 0.34 ± 0.04 keV and a 
best-fitting power-law photon index of 1.861q'[)3- 

Finally, we try to model the 0.3-10 keV X-ray spectrum 
in an entirely self consistent way. Therefore, we remove the 
two iron line species at 6.71 and 6.96 keV, represented by 
the two Gaussian components, and we add a second thermal 
component (Model 5b, Fig. [9l bottom panel), assuming in 
this way that these iron species originate from the emission 
spectrum of a hot diffuse gas. During the fit the normal- 
isations of the various model components are left to vary 
freely. The of the fit is similar to that of Model 5a, be- 
ing 900.70 for 767 d.o.f. (5.77 x 10"* NHP) yielding a low 
plasma temperature of fceTi ~ 0.36 ± 0.04, a high plasma 
temperature of fcBT2 ~ 8.S4'^\''[^ keV and a power-law in- 
dex of 1.86 ± 0.02. Note that this model is more reliable not 
only because it yields a marginally better NHP than Model 




1 2 3 

Observed energy (keV) 



Figure 9. EPIC spectral fits in the observed 0.3—10 keV energy 
range. (Top panel) Model 5a: Four Fe lines [Fe Ko (fixed), Fe xxv- 
r (fixed), Fe xxvi (fixed), Fe K/3 (fixed)]; a redshifted power-law; 
a thermal component. (Bottom panel) Model 5b: Four Fe lines 
[Fe Ka (fixed), Fe xxv-r (fixed), Fe xxvi (fixed), Fe K/3 (fixed)]; 
a redshifted power-law; two thermal components. The dotted and 
dashed lines correspond to the 2 and 4 per cent cross calibration 
uncertainties. 



5a, but it also provide us a physical explanation about the 
origin of the ionised iron line species. 

The best-fit results are plotted in the bottom panel of 
Fig.[9]having no significant residuals remaining in the 0.3-10 
keV energy band. The relatively poor value of is mainly 
due to the residuals, below 1.5 keV. The ratio plot, indicates 
that the pn, MOS 1 and MOS 2 residuals are consistent with 
each other at energies larger than 0.5 keV, and their ampli- 
tude is mostly smaller than 2 per cent in the whole energy 
range between 0.3 and 1.5 keV (and up to 4 per cent in a 
few cases; the dashed an dotted lines in this plot indicate 
the 2 and 4 per cent amplitudes, respectively, in the 0.3- 
1.5 keV band). On the other hand, the MOS 2 residuals do 



© 2002 RAS, MNRAS 000.[]M2l 



8 D. Emmanoulopoulos et al. 



Table 1. Best-fit spectral results for the EPIC science products. 



Model-parameter 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5a 


Model 5b 


[\o—0) U Kev 


o— lu Kev 


0— iu Ke V 


Q in \rr,\T 

vj— iu Ke V 


n Q in ^rci\T 
U.o— iU Kev 


n Q in 
U.o— lu Kev 


la) zpowerlw: Photon index, V 


1.66 ±0.02 


1.63 ±0.02 


1 66+°°^ 
^•""-0.02 


1 66+"°! 
^■""-0.02 


-, o-,+0.02 
^■*O-0.03 


1.86 ±0.02 


lb) zpowerlw: Normalisation^ 


2.84 ±0.01 


2.78 ±0.02 


2.87 ±0.02 


2.88 ±0.02 


3.27 ±0.02 


3.15 ±0.02 


2a) zGauss: Fe Ko line-energy (keV) 




6.42 ±0.01 


6.41 ±0.01 


6.41 ±0.01 


6.41 


6.41 


2b) zGauss: Normalisation^ 




1 09+0.16 
^■"^-0.20 


1 oc-i-o.ie 

l-8t>_0.19 


-1 70+0.15 


1 trc + O.lS 
^■'''^-0.16 


1.48 ±0.14 


3a) zGauss: Fe K/3 line-energy (keV) 




7.06 


7.06 


7.06 


7.06 


7.06 


3b) zGauss: Normalisation'^ 




0.14x2b* 


0.14x2b* 


0.14x2b* 


0.14x2b* 


0.14x2b* 


4a) zGauss: Fe xxv-r line-energy (keV) 






°- ''^-0.05 


6.71 ±0.02 


6.71 


— 


4b) zGauss: Normalisation'^ 






7 q-t-1.2 
' ■•^-2.6 






— 


5a) zGauss: Fe xxvi line-energy (keV) 










6.96 


— 


5b) zGauss: Normalisation'^ 










1 3+^-^ 





6a) mekal: Temperature, k^T (keV) 








0.34 ±0.04 


0.36 ±0.04 


6b) mekal: H density, njj (cm~^) 










1 


1 


6c) mekal: Metal abundance 










1 


1 


6d) mekal: Normalisation^ 










*-'^-1.30 


4.58l}:^8 


7a) mekal: Temperature, k^T (keV) 












8.84+};![* 


7b) mekal: H density, njj (cm~^) 












1 


7c) mekal: Metal abundance 














7d) mekal: Normalisation^ 












50.2+?l- 




404.57/387 


696.13/547 


617.05/544 


593.41/541 


896.44/764 


900.70/767 



Values without uncertainties correspond to fixed fitting parameters. 

The best-fitting normalisations correspond to the average of the two MOS (tied together) and pn. 
* 14 per cent of the Fe Ka line's normalisation value (parameter 2b). 
1 In units of 10~^ phot. keV'icm-^ s'l at 1 keV. 
^ In units of 10~^ phot. cm~^ s~l at the line energy. 
3 In units of J n^nndV, (lO"^^ cra'^). 



show discrepancies with the best-fit model [and the pn and 
MOSl) of the order of 6 per cent at energies below 0.5 keV, 
and in fact they contribute with a significant part to the 
resulting x^- Given the fact that the best model- fit residu- 
als do not exhibit well structured features (indicative of the 
presence of an extra spectral component that we have not 
taken into account) , we conclude that the low amplitude soft 
energy band residuals may be due to remaining calibration 
uncertainties of each detector, and also un-modelled cross- 
calibration uncertainties between the individual detectors as 
welS 

Finally, we also investigate whether the addition of a 
second, low energy thermal component could improve the 
quality of the model fit (previous X-ray spectra of the source 
have been modelled by two thermal components in the soft 
energy band; see the Discussion section). However, the addi- 
tion of this spectral component does not improve the quality 
of the fit, since it yields a of 899.12 for 764 dof (5.04x 10"* 
NHP). In fact, the best-fit temperature of the second com- 
ponent is consistent with zero (with a 90 per cent upper 
confidence limit of 9.7x10"^ keV). This was the case both 
when we leave the temperature of the first thermal compo- 
nent and the power-law spectral slope to vary freely during 
the model fit, and when we keep them frozen to their best-fit 
values of the previous model fit. 
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After estimating the best-fitting model parameters of 
Model 5b for the overall spectrum, we can compute the X- 
ray flux and X-ray luminosity of NGC7213. In the 0.5-10 
keV and 2-10 keV energy ranges the absorbed [unabsorbed] 
flux is (1.92(1.99] ± 0.01) x 10"" erg cm"^ s'l and 
(1.22[1.23] ±0.01) X 10~" erg cm"^ s"\ respectively. The 
corresponding X-ray luminosities for the same energy bands 
are (1.12[1.16] ±0.02) x 10''^ erg s'^ and (7.14[7.20] ±0.03) x 
10*1 erg s~i, respectively (the errors in the fluxes and the lu- 
minosities indicate the 68 per cent c.i. in order to be compa- 
rable with previous observations). Using the 2-10 keV lumi- 
nosi ty, we can estimate t he expected EW of Fe Ka using eq. 
2 in lBianchi et al ] (I2OO9'). This yields a value of around 124 
eV being very consistent with our estimate, being around 
126 eV (Section 112111. 



5 RGS X-RAY SPECTRAL ANALYSIS 

We first model the 0.35-1.8 keV RGS spectra of NGC7213 
independently on the results of the EPIC spectral fitting, 
in the overlapping energy ranges (i.e. above 0.4 keV). The 
best-fit parameters for the following X-ray spectral models 
are given in Table [2] 

Our initial continuum model includes a power-law at- 
tenuated by neutral absorption, A'h (model A). During the 
fit, the power-law spectral indices between the two RGSs are 
linked together in contrast to the power-law normalisations 
which are left to vary independently between the two RGSs, 



|http : //yam ■ esac ■ esa ■ Int/ext 6rnal/xmm_user_support/usersgroup/2(t tgadaQ)4litdfar. flciaslble uncertainties in the cross-calibration 
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between the two spectrometers. The fit is statistically ac- 
ceptable with a of 298.59 for 269 d.o.f. (0.10 NHP) yield- 
ing a best-fitting power-law photon index of 1.78 ± 0.03. 
However, the ratio plot for both RGSs (Fig. 1101 top panel) 
show broad systematic deviations from a constant, as well 
as two line-like features at around 19 A (0.653 keV) and 22 
A (0.564 keV, only for RGS 1), in the observer's frame. 

In order to cure both the 10-20 A and 25-30 A broad 
excesses in the residuals, and at the same time to model the 
two line-like emission features, we add to the intrinsic nu- 
clear continuum the mekal-model, parametrised by its tem- 
perature and normalisation (emission- measure of the gas). 
Similarly as before both the Hydrogen volume density and 
the chemical abundance of metals are set to unity and to 
that of Solar abundance, respectively. During the fit, we 
link together the temperature of the plasma for the two 
RGS spectra, while leaving their normalisations to vary in- 
dependently between them. As for the power-law, we fix the 
photon index to the best-fitting value derived by Model 5b, 
1.86 (Sect. f43)) . since the EPIC data trace better the source 
continuum due to their high sensitivity and energy coverage. 
This RGS best-fitting model (model B) describes the data 
significantly better than the simple power-law best-fitting 
model having a of 261.85 for 267 d.o.f. (0.58 NHP) and 
yielding a best-fitting temperature of 0.39 ± 0.07 keV. The 
best-fitting model does not leave any significant residuals, 
associating at the same time the emission lines from the ther- 
mal plasma, at 19 and 22 A, with a prominent OvillLya 
and OvilKa-f, respectively (Fig. 1101 bottom panel). The 
best-fitting normalisations of the mekal components for the 
two RGSs are consistent between each other, within their 68 
per cent confidence limits, and both of them are also fully 
consistent with the best-fitting EPIC values. In addition the 
best-fit temperature is entirely consistent with that of Model 
5b using the EPIC data. From a statistical point of view the 
addition of the t hermal component can be jus tified by per- 
forming an f-test (|Bevington fc Robinson|[l992l ) between the 
best-fitting model A and model B yielding a probability of 
1.44 X 10~® (for an f-statistic value of 14.143) indicating that 
the addition of the thermal component is valid for 1 per cent 
significance level. 

The above best-fitting model describes adequately the 
RGS data of NGC 7213. However, due to the limited RGS 
spectral range coverage, 0.35-1.8 keV, and to the lack of 
sufficient sensitivity at the OvilKa and NelXKa triplel[f| 
wavelengths, RGSs data alone can neither constrain strongly 
the existence of plasma with temperatures significantly 
higher or lower than that of our best-fitting RGS tempera- 
ture at around 0.3 keV. 



^ The RGS 1 is missing the spectral region where Ne IX linos 
are usually located (between 11-14 A), and the RGS 2 is miss- 
ing the region around the Ovil lines (between 20-24 A). Both 
instruments are heavily affected from bad columns. In order to 
quantify this, we count the number of independent spectral ele- 
ments {5X = 60 mA) for each RGS affected by bad columns, and 
we divided it by the total number of independent elements in the 
analysed bands yielding a fraction of 14 per cent of bad columns 
for the RGS 1 and 12 per cent for RGS 2. 
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Figure 10. RGS spectral fits in the observed 0.4—1.6 keV energy 
range. The vertical labels indicate the energy in the source's rest- 
frame. (Top panel) Model A: A redshifted power-law. (Bottom 
panel) Model B: A thermal component; a redshifted power-law 
(fixed). 



6 SUMMARY AND DISCUSSION 

We have analysed the longest XMM-Newton observation of 
NGC 7213, obtained during November 2009. We have per- 
formed hardness ratio analysis, on the EPIC light curves in 
different energy bands, yielding no evidence of X-ray spec- 
tral variability within the given observation. This enabled 
us to perform a detailed X-ray spectral analysis on both 
the EPIC and the RGS data products using the total X-ray 
spectra, in the 0.3-10 keV and 0.35-1.8 keV energy bands, 
respectively. 

The main results from our X-ray spectral analysis can 
be summarised as follows. 
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Table 2. Best-fit spectral results for the RGS science products in the [(0.35 - 0.9) U (1.2 - 1.8)] keV. 



Model-parameter Model A Model B 

la) zpowerlw: Photon index, V 
lb) zpowerlw: Normalisation^ 
2a) mekal: Temperature, k^T (keV) 
2b) mekal: H density, nj] (cm~^) 
2c) mekal: Metal abundance 
2d) mekal: Normalisation^ 

X^/d.o.f. 289.59/269 261.85/267 



Values without uncertainties correspond to fixed fitting parameters. 

The best-fitting normalisations correspond to the average of the two RGSs. 

* Individually the normalisations for RGS 1 and RGS 2 are 7.47 ± 2.12 and 4.93^^ g|^ respectively (in the corresponding units). 
^ In units of 10~^ phot. keV~^cm~^ s~^ at 1 keV. 
2 In units of J n^n^dV , (10^^ cm-3). 



1.78 ± 0.03 
3.10 ±0.05 



1.86 
2 07+0.07 

-0.05 

0.39 ±0.07 
1 
1 



6.21 



-1-2.21 , 

-2.07 



• The EPIC X-ray spectrum: 

- above 3 keV exhibits three narrow emission iron lines 
at 6.41 ± 0.01, 6.71 ± 0.02 and 6.96t|3;!il^ being consistent 
with the neutral Fe Ka (EW: 126+2? eV), the Fe xxv 
resonance (EW: 58^20 eV) and the ls-2p doublet Fe xxvi 
Ka (EW: 23l^^ eV). 

- in the full band (0.3-10 keV) can be best-fitted by 
an underlying power-law model with F = 1.86 ±0.02, and 
two thermal models, (mekal), of temperatures equal to 
fcBTi = 0.36 ± 0.04 keV and kBT2 = 8.84+^1* keV (with 
the highest temperature component effectively accounting 
for the ionised iron emission lines) 

• The RGS X-ray spectrum: 

- The results from the RGS X-ray spectrum analysis in 
the 0.35-1.8 keV band are fully consistent with the best- 
fitting model of the broad band EPIC spectrum. In par- 
ticular, the RGS spectrum can be well fitted by a power- 
law component and a mekal component of temperature 
fceT = 0.39 ± 0.07 keV. 

The best-fitting photon spectral index that we get from 
the EPIC data together with the fiux estimates (Model 5b, 
Sect l4.3p . are fully consistent, within the 68 per cent confi- 
de nce limits, with the linear anticorrelation relation derived 
bv lEmmanoulopoulos et al.l l|2012l ). The latter is deduced us- 
ing all the archival RXTE data, showing unambiguously that 
NGC 7213 behaves similarly to Galactic 'hard state' sources 
exhibiting a 'harder when brighter' X-ray behaviour. 

The EW value of the neutral Fe Ka is significantly 
larger than the one derived from e arlier XMM-Newton 
observations (around 82 ± 10 eV; iBianchi et al.l l2003l : 
[star ling ct al. 20 0j) and S uzaku observations (around 83 ± 
10 eV: lLobban et al.ll2010l ) (note that the error bars indicate 
the 90 per cent confidence limits of the corresponding EWs). 
However, it is absolutely in accordance with the one derive d 
from the Chandra-RETG {UOtfo eV: iBianchi et al.ll2008l ). 
In the case of X-ray reflection from Compton-thick mate- 
rial in a plane-like configuration, the EW of the narrow iron 
emission line should be tightly connected with the refiection 
strength R. In fact, one considers a linear correlation be- 
tween the two parameters, b eing equal t o unity whe n the Fe 
Kg li ne's EW is 120 eV (e.g. lMatt et al.|[l991 : Bianchi et all 
l2009l ). then the EW value of the neutral Fe Ka line that we 



measure implies a refiection strength of i? ~ 1, correspond- 
ing to an average solid angle, subtended by the Compton- 
thick refiector of 27r sr. This is not consistent with our results 
disfavouring the Compton-thick refiector origin of neutral Fe 
Ka . 

This result favours scenarios in which the neutral Fe 
Ka line origi nates from a C ompton-thin material, such as 
th e BLR (e.g. IBianchi et al.| [2008). Following the procedure 
of lYaqoob et all (l200lD ~ assuming a spherically symmetric 
cloud distribution and neutral iron, the expected EW for 
th e BLR in NGC 7213 can be written (in a similar fashion 
to IBianchi et al.l l|2008l )) as 



EWf 



= 31.7 I 



(0.35) ( 



1023 c 



eV 



(1) 



where /c is the covering factor and A^h the hydrogen col- 
umn density of each cloud. For this relation we assume an 
iron abundance relative t o hydrogen of Apc = 4.68 x 10~^ 
jAnders fc Grevess3ll989l ) and the fiuorescence yield is set 
to uj]^ = 0.342 for neutral iron (iBambvnek et all 1 19721 ). As- 
summg fc = 0.35 (|Goad fc Koratkarll 19981 ) a column density 
of A'^h = 4 X 10^^ cm~2 is required to reproduce the Fe Ka 's 
EW. Thes e values of fc and A'h are commonly expected from 
the BLR (|NetzeJI 19901 ) ■ making the latter indeed a possible 
candidate for the origin of the neutral Fe Ka emission line 
in NGC 7213. 

The X-ray spectrum of NGC 7213 is known to ex- 
hibit highly ionised iron emission lines between 6.5-7 keV 
tei anchi fc MattI 120 02: Bian chi etahl '2005: St arling et"al] 
120051 : IBianchi et al.i ,2008 ; .Lobban et al.. ,2010. ). From our 
EPIC data, we find that the dominant component of the 
He-like Fe xxv triplet is the resonant one. This is in agree- 
ment with both older XMM-Ne wton data ([Star ling ct aP 
I2OO5I ) and Cfeanrfra-HETC data llBianchi et aiT[2008l ) but 
not with the Suzaku data ([Lobban et al.ll2010l ). the latter 
associating it with the forbidden transition. The resonant 
nature of Fe xxv line indicates th at it is produ ced within 
a gas which is coUisionally excited ([Porauet fc D ubau 200_3; 
iBautista fc KallmanI [19991 ). 

Finally, both the EPIC and the RGS data suggest the 
presence of a multi-temperature hot gas in the nuclear region 
of NGC 7213. Our broad band EPIC resuffs show that the 
source's 'soft excess' component, at energies below 2 keV, 
can be well fitted by one thermal component with best- 
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fitting temperature of around 0.36 keV, something wliicli 
is in accordanc e witli the RGS dat a. Previous results of 
NGC 7213 from iLobban et aP (|2010l ) consider two thermal 
emission components to best-fit the l ow-energy part of the 
X-ray spectrum at 0.24 and 0.86 keV. IStarling et"all (120051 ) 
found evidence for a single low energy thermal component 
at 0.18 keV and by fitting a second one, at 0.56 keV, the fit 
was not significantly improved. Both the pn and MOS spec- 
tra from the long XMM-Newton observation that we analyse 
in this work are not consistent with the presence of a second 
thermal component. Interestingly, our best-fit temperature 
of the single low-energy t hermal component is intermediate 
between that reported by iLobban et al.1 (|2010l ). This could 
perhaps explain the difference between ours and the previ- 
ously reported results. Finally, the thermal component with 
the higher temperature, 8.84 keV, can effectively account for 
the previously mentioned ionised iron emission lines, that we 
detect between 6.5-7 keV, being indicative for a gas in col- 

lisional ionisation equilibrium. 

The Chandra X-ray image of the source l|Bianchi et al.l 

I2OO8I ) shows only an unresolved nuclear source (even at low 
energies) within a few arcseconds around the galactic nu- 
cleus. Consequently, if a multi-temperature gas does exist 
in this source, it has to be confined in a small region, close 
to the AGN. Perhaps then, intense starburst superwind ac- 
tivity in the nuclear region of this source could account of 
this multi-temperature medium, which may even be outflow- 
ing at velocities of the order of 1000 kr n s~^, as suggested b y 
the iron line's blueshifts measured bv lBianchi et al.l (|2008t ). 
In the framework of superwinds, we expect the hard X- 
ray emission from the galaxy driving the superwind com- 
ing from the hot and tenuous supernova/wind-heated gas 
inside the starburst (i.e. inside the son ic radius of the wind) 
l|Heckman et al.l ll993: Strickland 2001). In this environment, 
characteristic temperatures between 6-9 keV can arise yield- 
ing in a self-consistent way the observed ionised iron species. 
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